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Abstract 

Ambiguities in the definition of angular momentum of a quantum-mechanical 
particle in the presence of a magnetic vortex are reviewed. We show that the 
long-standing problem of the adequate definition is resolved in the framework of 
the second-quantized theory at nonzero temperature. Planar relativistic Fermi 
gas in the background of a point-like magnetic vortex with arbitrary flux is 
considered, and we find thermal averages, quadratic fluctuations, and correlations 
of all observables, including angular momentum, in this system. The kinetic 
definition of angular momentum is picked out unambiguously by the requirement 
of plausible behaviour for the angular momentum fluctuation and its correlation 
with fermion number. 
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1 Introduction 



Whereas in classical theory a comprehensive description of electromagnetism is given 
in terms of the electromagnetic field strength acting locally and directly on charged 
matter, this is not the case in quantum theory. Since the quantum-mechanical wave 
equation for a charged particle involves the electromagnetic vector potential, the state 
of charged matter can be influenced by electromagnetic effects even in the situation 
when the space-time region of nonvanishing field strength is not accessible to charged 
matter P3I2|- The latter was demonstrated for the case of quantum-mechanical scat- 
tering on a tube of the magnetic force lines, which was impenetrable for scattered 
particles: the differential scattering cross section was found to be a periodic function 
of the total magnetic flux confined in the tube (2jj. This phenomenon with all its var- 
ious generalizations (see, e.g., Refs.jHl IH E]) has no classical analogues and is denoted 
usually as the Bohm-Aharonov effect. 

Fractionalization of quantum numbers is a well-known feature of various field theory 
models that describe interaction between fermions and topological solitons 13 El 13 
Uni im CSj . It is usually convenient to interpret solitons as background fields and to 
quantize fermion fields in such backgrounds. In this respect it seems to be of interest 
to consider a case when the region where a background field is nonvanishing does not 
overlap with the region where a fermion field is quantized. Namely, the quantized 
fermion matter does not penetrate a tube of the magnetic force lines; for brevity, the 
impenetrable magnetic tube will be called the magnetic vortex in the following. One 
inquires about the dependence of the properties of this second-quantized system on the 
vortex flux and the boundary condition on the edge of the vortex; if a transverse size 
of the vortex is neglected, then a parameter of the boundary condition exhibits itself 
as a parameter of a self-adjoint extension of the one-particle hamiltonian operator. 

A study of this second-quantized system started two decades ago [inillll. It was 
shown for a particular choice of the boundary condition that fermion number ^3] , mag- 
netic flux and angular momentum ^Tj are induced in the vacuum of quantized 
massive fermions. The induced vacuum quantum numbers under the most general 
set of boundary conditions were obtained in Refs. [THl IT^ 1201 12I|- All results, except 
angular momentum, are periodic in the value of the vortex flux, and such a type of 
behaviour is consistent with the above mentioned for the quantum-mechanical Bohm- 
Aharonov effect. As to angular momentum, the situation is not so far determinative, 
which is due to two alternatives for its definition. In quantum field theory, the angu- 
lar momentum operator, as well as operators of all other observables, is obtained by 
sandwiching the appropriate quantum-mechanical (one-particle) operator between the 
fermion field operators (see Section 4 below). In quantum mechanics, there are two 
possibilities to define angular momentum of a fermion particle in the background of the 
Bohm-Aharonov magnetic vortex: the canonical one that is quantized with half-integer 
eigenvalues. 
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and the kinetic one that is quantized otherwise, 

J = n + ^ - e$, (2) 

where e is tlie electric charge and $ is the vortex flux in units of 27r; strictly speaking, $ 
is measured in units of 2'irhc, but we use conventional units h = c = 1. The motivation 
in favour of the canonical definition lies in that it yields angular momentum which 
is conserved when the vortex flux is varied in time [21 1221 1221 121] ■ On the other 
hand, the kinetic definition proves to be rather rewarding, since it leads to such prolific 
concepts as anyons and fractional statistics [2312^1 which are of fundamental value and 
possess fascinating phenomenological applications to some condensed matter systems 
j27| l28j l29t 180] . It should be noted that quantum-mechanical consequences of the 
kinetic definition are periodic in the value of e$ with the period equal to 1, and this 
can be regarded as one more of manifestations of the Bohm-Aharonov effect. If we 
turn to quantum field theory and consider angular momentum which is induced in the 
vacuum around a magnetic vortex, then the kinetic definition yields the quantity which 
depends on the fractional part of e$, while the canonical definition yields the quantity 
which depends on both the fractional and integer parts of e$ [13 |2I] • 

In the present paper we shall study the effect of nonzero temperature on the induced 
angular momentum around a magnetic vortex, and consider its thermal average and 
quadratic fluctuation, as well as correlations with other observables. The previous 
analysis in Refs. |3H lH^ is augmented by taking account of the bulk volume contribution 
to thermal characteristics. In general, the expressions for thermal characteristics consist 
of two pieces: the one corresponding to the ideal gas contribution depends essentially 
on a size of the system, increasing by power law as the size increases, and the other one 
corresponding to the correction due to interaction with a magnetic vortex is finite in the 
limit of the infinite size; although, sometimes, as are the cases of the averages of fermion 
number and induced magnetic flux and their correlation, the ideal gas contribution is 
vanishing. We shall show that the canonical definition of angular momentum yields 
rather unconvincing results for the angular momentum fluctuation and its correlation 
with fermion number: these quantities contain pieces that are proportional to both 
powers of the system size and powers of the whole vortex flux, and such a mismatch 
of different contributions is hardly tolerable. On the contrary, the kinetic deflnition 
of angular momentum yields plausible results, and this allows one to conclude that 
quantum fleld theory at nonzero temperature favours for certain the latter deflnition. 

In the next section, we review the deflnition of angular momentum for a quantum- 
mechanical charged particle placed into a classical static rotationally invariant magnetic 
fleld. In Section 3 we consider a situation when the particle does not penetrate the 
magnetic fleld region and discuss some manifestations of non-simply-connectedness 
of the accessible to the particle region. A general formalism for the treatment of 
thermal characteristics of observables in the framework of quantum fleld theory at 
flnite temperature is introduced in Section 4. Section 5 deals with planar relativistic 
Fermi gas in the background of a magnetic vortex, and we determine relevant one- 



3 



particle spectral densities. Using formalism of Section 4 and results of Section 5, we 
obtain thermal characteristics of all observables in this system in Section 6. We discuss 
our results in Section 7. Some details in the derivation of the results are outlined in 
Appendices A-D. 

2 Angular momentum of a particle in a magnetic 
field 

Angular momentum of a free particle consists of the orbital and spin parts 

Jo = -ix X 9 + i(T. (3) 

If a charged particle is posited into an external electromagnetic field, then a contri- 
bution of the latter has to be added to Eq.(3). It can be shown (see Appendix A, 
Eq.(A.18)) that this contribution in its turn consists of the purely field part 

Mem = y'rfV(x'x(ExB)), (4) 

and the particle-field interaction part 

M/ = -ydVx'x A/ = -exx A(x), (5) 

where 

/ = e5=^(x'-x) (6) 

is the charge density of a particle located at point x. Adding Eq.(5) to Eq.(3) results 
in the emergence of the covariant derivative, d ^ d — ieA(x), in the orbital part. 

Let us dwell now on the purely field contribution to the angular momentum, Eq.(4). 
Substituting B = 9 x Ay (where Aj is a certain, fixed, vector potential) into Eq.(4) 
and integrating it by parts, we get 

+ j d?x'e'^'^x'\AJd'^E'' + E'^d'^A'} + A^id'^E"^ - d'^E"")]. (7) 

Note that Mem is the improved angular momentum which differs from the canonical 
angular momentum Mem (see Eqs.(A.ll) and (A. 12)), 

M'em - M'em = / d^x' d[s'^'^x'^x'^'^^ = - j d^x' ^is'^'^x'^ A'^E^ (8) 

where we have used Eq.(A.16) for x^^^- We consider a static rotationally invariant 
magnetic field, (x x ^) • B = 0, and choose Af satisfying conditions 

x-A/ = 0, 9-A/ = 0. (9) 
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As to the electric field, it is generated by static charge density j° (6), 

d'-^=f, E = f^^, (10) 

477 |x' — Xl'^ 

and is curlless, 9' x E = 0. Then only the first term in square brackets in the second 
integral in the right hand side of Eq.(7) survives, yielding upon integration a contribu- 
tion which is equal to the contribution of particle-field interaction, Eq.(5), taken with 
opposite sign and at A = A/. As to the first integral in the right hand side of Eq.(7), 
it is transformed to an integral over a closed surface, as described in Appendix A, see 
Eqs.(A.4) and (A. 13). Thus the behaviour of terms in brackets in this integral is cru- 
cial: the first two terms decrease as \^'\^^ at large |x'| yielding vanishing contribution 
upon integration, whereas the last one decreases as |x'|~^ at large |x'|. This yields upon 
integration finite contribution which is equal to the total magnetic flux times e(27r)~^ 
with the opposite sign [21]. Consequently, one gets 

MsM = exx A^(x)-;^ /rfVB, (11) 

Zn J 

where the integral is over a plane which is orthogonal to B, and 

MEAf = exx A/(x). (12) 
Summing Eqs.(3), (5) and (12), one gets 

J = -ix X [d- ?eA(x) + ieAf{x.)] + -cr. (13) 
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Taking into account relation 



X X A.(x) = — / rfVB, (14) 

D 

where the integration is over region D defined by {x'y + {x'^f < {x^f + {x^f, one 
can verify that J (13) commutes with the one-particle hamiltonian operator; thus J 
(13) can be regarded as the operator of total angular momentum of the system, which 
is conserved as a consequence of rotational invariance. 

If the improved quantity, Eq.(ll), is taken instead of the canonical one, Eq.(12), 
then, summing Eqs.(3), (5) and (11), one gets 

J = -i X X [9 - ieA(x) + ieA/(x)] - — / ciV B + -cr. (15) 

27r J 2 

Since the difference between operators J (13) and J (15) is x-independent and diagonal 
in spin indices, the latter operator commutes with the hamiltonian as well. 

Although both the canonical and improved definitions of angular momentum of 
electromagnetic field are compatible with the conservation of total angular momentum. 
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the motivation in favour of the first one may be seen in the fact that we start initially 
from the canonical angular momentum. Then, at intermediate stage, the improved 
angular momentum is introduced as a convenient and well-defined quantity which is 
constructed directly from the gauge invariant and symmetric energy-momentum tensor, 
see Appendix A. This quantity plays a somewhat auxiliary role allowing one to express 
the canonical angular momentum through the explicitly gauge invariant entities. Note 
that angular momentum J (13) in the gauge A = Af coincides with angular momentum 
of a free particle, Eq.(3). Thus, it is evident that spectra of operators Jo and J coincide, 
and this may serve as an additional argument in favour of the canonical definition of 
angular momentum. 

Note also that the difference between J (15) and J (13) is equal to the return 
magnetic flux (times e(27r)~^), see Eq.(8), which completes a flux loop at extremely 
large distances outside the position of the charged particle. If the particle is allowed to 
go beyond the region of the return flux, then, as it enters this region from outside and 
crosses it, one arrives at expression (13) for total angular momentum inside the region 
of the return flux. 

In the case of a relativistic massive spinor particle, considered in the present paper, 
the spin matrix is cr = {2i)~^a x ct and the hamiltonian is 

H = -ia-[d -ieA{^)] + -f°m, (16) 

where a = 7*^7, and 7°, 7 are the Dirac matrices. 

3 Bohm— Aharonov configuration of the 
magnetic field 

Of special physical interest, as it was noted in Introduction, is the case when the 
region of the magnetic flux lines is not accessible to a charged particle [21 • Namely, 
one considers a rotationally invariant magnetic fleld conflguration in the form of a long 
solenoid with extremely distant return flux, and the charged particle is posited in the 
fleld-free region between the solenoid and the return flux. From now on we shall call 
the region of solenoid as the inner region and the fleld-free region as the outer region, 
implying that the return flux is all pushed out to inflnity. 

In the outer region, improved quantity (4) is obviously vanishing, NLem = 0, 
whereas canonical quantity is given by the same expression as Eq.(12) but is 

x-independent, since in this case it depends on the whole flux through the inner re- 
gion, X X A/-(x) = 2^ / d'^x''B. Eq.(13) takes form 

J = -i X X [d - zeA(x)] + ^ J d^x'B + x a, (17) 

whereas Eq.(15) takes form of the kinetic (or mechanical) angular momentum of the 
particle, 

J = -i X X [d - ■ieA(x)] + —cx x ex. (18) 
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Note that the outer region is not simply connected, and one may feel free not to require 
the coincidence of the spectrum of total angular momentum with that of Jq (3). Also, 
vanishing of the improved angular momentum of electromagnetic field in this case may 
be regarded as an argument in favour of choosing J (18) in the capacity of total angular 
momentum of the system. Whether J (17) or J (18) gives the physically meaningful 
angular momentum, was disputed in the literature, see, e.g., Refs.|221 ESI SUES]- In the 
present paper we shall show that namely the kinetic definition of angular momentum 
provides the physically reasonable behaviour of the thermal quadratic fluctuation and 
correlation in the second-quantized theory at nonzero temperature. 

But before going to the second-quantized theory, let us discuss some manifestations 
of non-simply-connectedness of the outer region. Let magnetic field B be directed along 
the x'^-axis in the inner region, and the outer region correspond to (x^)^ + (x^)^ > 
with Tc being the transverse size of the inner region. Then vector potential Aj takes 
form 



2 1 

3y 3j 



A\ = -$ - Al = $ - Al = fl9) 

where 

$ = — / d^x'B^ (20) 



27r 

is the total magnetic flux in units of 27r, and Eq.(18) in gauge A = A/ takes form 

J3 = -i(x^d2 - x'^di) - e$ + ^a^a^ (21) 

^ t 

Operator (21) acts on functions which are single-valued in the outer region: 

(r, x^l = (r, (^ + 27r, x^l, (22) 

where r and (p are the polar coordinates in the x^x^-plane. Since the outer region is not 
simply connected, one can consider functions satisfying much more general condition 

(r, x^\' = e'^^'^ir, ip + 27r, x^\', (23) 

where S is a continuous real parameter. Functions satisfying conditions (22) and (23) 
are related by gauge transformation 

(r, (f, x^\' = e"*^'^(r, (f, x^\, 

A' = Af- e'^ E d<^. (24) 

Neglecting the transverse size of the inner region (rc 0), and presenting the angular 
variable with range < (/9 < 27r as 

ip = arctan ( ^ j + n [e{-x^) + 2^(x^) ^(-x^)] , (25) 
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where 9{u) — | g' < g |' S®^^ immediately 

and, consequently, 
1/1 



A'' = -($-e-^S) 



^'^ = - ^"-) (^T)^T(^ + ^^e-^mx')6{x% (27) 
A" = 0. 

Although the curl of the vector potential remains invariant under gauge transformation 
(24), 

^3 = diA} - d2A) = ^1^'' - d2A'^ = 27r$(5(a;i)(5(a;2), (28) 
angular momentum (21) is changed to 

J'3 = -i{x^d2 - x^di) - e$ + S + V^ (29) 

2z 

but its spectrum remains unchanged, because operator (29) acts on functions satisfying 
condition (23). In particular, choosing S = e$ one gets the gauge with the vector 
potential eliminated everywhere excepting the semiaxis of positive x^: 

A'' = 27r$^(x^)5(x2), A'^^0. (30) 

Kinetic angular momentum (18) in this gauge coincides with angular momentum of a 
free particle, Eq.(3), but their spectra differ: the operator of the latter acts on functions 
which are defined everywhere and, thus, single- valued (22), while the operator of the 
former acts on functions which are defined in space with a cut along the positive x^ 
semiaxis and, thus, satisfying boundary condition (23) with S = e$ on the sides of 
the cut. The cut can be continuously deformed, and, in general, its projection on the 
x-'^x^-plane can be a curved line starting from the origin and going to infinity. Eq.(30) 
in the case of such cut takes form 



A'' = -27r$ I e^"'S^{yi - x')dx'^\ 



where coordinate vectors x and x' lie in the x^x^-plane. 

To conclude this section, we note that, obviously, canonical angular momentum 
(17) in gauge A = A/ takes form 

J3 = -i{x^d2 - x^di) + ^aV^ (31) 
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and it is changed to 

J'3 ^ -i(x^d2 - x'^di) +E+ ^a^a^ (32) 

under gauge transformation (24). The spectrum of canonical angular momentum is 
given by Eq.(l), while the spectrum of kinetic angular momentum is given by Eq.(2). 



4 Second-quantized theory at finite 
temperature 

The operator of the second-quantized fermion field in a static background can be pre- 
sented in the form 

^(x,t) = ^ e-^^*(x|A)a, + ^ e-^^*(x|A)6+ , (33) 

{Ex>0) (Ea<0) 

where aj" and ax {b^ and bx) are the fermion (antifermion) creation and destruction 
operators satisfying anticommutation relations, 



ax, ay 



= bx,bt,\^ = {X\X'), (34) 

and (x|A) is the solution to the stationary Dirac equation, 

H{^\X) = Ex{^\X) , (35) 

H is the Dirac (one-particle) hamiltonian, A is the set of parameters (quantum numbers) 

specifying a one-particle state, and Ex is the energy of the state; symbol ^ means the 

summation over discrete and the integration (with a certain measure) over continuous 
values of A. Ground state |vac) of the second-quantized theory is defined as 

a^lvac) = bx\v8.c) — 0. (36) 

Let J be an operator commuting with the hamiltonian in the first-quantized theory, 

[J,H]_ = Q. (37) 

In the case of unbounded operators, commutation of their resolvents is implied, or, to 
be more specific, it is sufficient to require that operators H and J have a common set 
of eigenf unctions, i.e. relation 

J(x|A)=ja(x|A) (38) 

holds as well as Eq.(35). Eigenfunctions (x|A) satisfy the conditions of completeness 
and orthonormality; in general, normalization to a delta function is implied. Thus, 
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in the second-quantized theory, the operators of the dynamical variables (physical 
observables) corresponding to H and J can be diagonalized: 



U=]^j d'^x [vl/+(x,t), /7^(x,t)]_ = ^E, 



and 



M = 



(fx 



vl/+(x,t), Jvl/(x,t) 



«A «A - &A ^A - -Sgn(£;A) 



«A - &A ^A - -Sgn(EA) 



(39) 



(40) 



d is the space dimension, and sgn(M) = 9{u) — 9{—u) is the sign function. 

The thermal average of the observable corresponding to operator (40) is conven- 
tionally defined as (see, e.g., R,ef.|HH]) 



M{T) = (M), ^ : _^,y , (5 = {kBT)-\ 



Spe 



-f3U 



(41) 



where T is the equilibrium temperature, ks is the Boltzmann constant, and Sp is 
the trace or the sum over the expectation values in the Fock state basis created by 
operators in Eq.(34). Also one defines the thermal quadratic fluctuation 



A(r;M,M) = (M2)^-((M)^)2. 



(42) 



Quantities (41) and (42) can be expressed through the derivatives of the appropriate 
thermodynamic potential: 



M{T) 



11=0 



A{T;M,M) 



1 d^n{(3,fi) 
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^l=a 



(43) 



where /i is the appropriate chemical potential, and 



fi(/?,/i) = --ln5pexp[-/5(f/-/iM)]. (44) 

We show in Appendix B that thermodynamic potential (44) is presented as 

f)(/3,/i) = --iTrlncosh[l/?(i7-/iJ)], (45) 

where Tr is the trace of an integro-differential operator in the functional space: Tr . . . = 
/ (i'^xtr(x| . . . |x); tr denotes the trace over spinor indices only. Then, using Eq.(43), 
one can express average (41) and fluctuation (42) through functional traces of operators 
in the first-quantized theory: 



M(T) = --TrJtanh(-/5i7) 



(46) 
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and 

A(T; M, M) = ^Tr J2 sech\-PH). (47) 
Eqs.(46) and (47) are transformed into integrals over the energy spectrum, 

1 °° 1 
M{T) = -2 / dETj{E) t^nh{-PE) (48) 



and 

where 
and 



1 °° 1 
A{T;M,M) = - I dETAE)sech'{-PE), (49) 



Tj{E) ^TrJS{H - E) (50) 



rj2{E) ^TrJ'^5{H - E) (51) 

are the appropriate spectral densities. 

Note that in the case ol J — I, where / is the unit matrix in the space of Dirac 
matrices, the corresponding operator in the second-quantized theory is the operator of 
fermion number; and Vt are the conventional chemical and thermodynamic potentials 
in this case. In the d = 1 case fermion number is the only observable which is conserved 
in addition to energy. In more than one dimensions there are more conserved observ- 
ables. In particular, in the d — 2 case, in addition to energy and fermion number, also 
total angular momentum is conserved when the system is rotationally invariant, as it 
was discussed in two previous sections. 

If an observable is not conserved, then its operator in the first-quantized theory does 
not commute with hamiltonian, [T, H] ^ 0, and its operator in the second-quantized 
theory, 

d = ^|d'^x[^+(x, t), Tvl/(x,t)]_, (52) 

is not diagonalizablc. Similarly to Eq.(48), the thermal average of the nonconserved 
observable can be presented as 



0{T) - -2 / dETr{E)i8.nh{-PE), (53) 
—00 

where 

Tr{E) = T:iT 5{H - E). (54) 

Also one can consider the thermal correlation of the conserved and nonconserved ob- 
servables, 

A(T; 6, M) = (6 M)0 - {d)p{M)p, (55) 
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which, similarly to Eq.(49), can be presented as 



1 °° 1 
A(T; d,M) = - I dETrjiE)sech\-PE), (56) 

— oo 

where 

Trj{E) = TtTJ6{H - E). (57) 

It should be noted that relations (45), (48), (49), (53), and (56) are somewhat 
formal, and the proper treatment of the normal ordering of operator product in the 
second-quantized theory is required. In the absence of interaction, the operators are 
normal ordered (see, e.g., Ref.|34j), i.e. the c-number pieces in Eqs.(39) and (40) 
are dropped. Consequently, in the presence of interaction with external fields, the 
c-number pieces in Eqs.(39) and (40) are renormalized by subtracting these dropped 
pieces. Correspondingly, thermodynamic potential (45) is decomposed as 

fi(/3,/i)=fiW(/5,/i) + fi«(A/i), (58) 

where 

fi(°)(/?,/i) = -iTrln{l + exp[-/3(|i7o| -Msgn(//o))]} (59) 
is the thermodynamic potential in the free field case, and 

■1 



f](i)(/3,/i) = -i jxrlncosh 



\(3{H-^,J) 



— Tr In cosh 



(60) 



is the addition which is due to interaction with external fields; here the operators in 
the first-quantized theory without interaction are denoted by Hq and Jq. Note that 
term — |Tr [\Ho\ — /iJoSgn(ifo)] is dropped, which corresponds to the normal ordering 
at zero temperature in the free field case. 
As a consequence of Eqs.(58)-(60) one gets 

M(T) = M(°) (T) + M(^) (T) (61) 

and 

A(T; M, M) = A(°)(T; M, M) + A(^)(T; M, M), (62) 

where 

oo , , 

M(°)(T)= /rfErf(E)J^, (63) 

— oo 

M«(T) = -- I dE4'\E)te^nh{-PE), (64) 

— oo 

Tf\E) = Tr Jo 6{Ho - E), Tf\E) = Ti J 6{H - E) - Tr Jq 5{Hq - E), (65) 
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and 



1 °° 1 
A(°)(r;M,M) = - / dETf^{E)sech\-pE), 



(66) 



A(^)(r; M,M) 



1 



(67) 



rS?(E) = TrJ^ 5(i7o - E), t)',' {E) = TrJ' 5{H-E) - Tr S{Ho - E). 



-(1) 



(68) 



Using relation 6{H — E) = [{H — E — iO)^^ — {H — E + iO)^^], one can trans- 
form integrals over the real energy spectrum in Eqs.(64) and (67) into integrals over a 
contour on the complex energy plane, thus yielding a representation of thermal char- 
acteristics through the renormalized resolvent traces, 



c 



and 



AW(r;M,M) = - / — : Trj2(//-a;) 



du 



c 



2m L 



^2 



(69) 



(70) 



where C is the contour consisting of two collinear straight lines, (— oo + iO, +oo + iO) 
and (-I-00 — iO, — oo — iO), in the complex a;-plane, and the renormalized resolvent traces 
are 

\VyJ{H - ^y\^^ = Tr J{H - u)-^ - TrJo{Ho - o;)"^ (71) 

and 



,-1 



(72) 



To conclude this section, we note that thermal average (53) and correlation (56) 
can be treated in a similar way. 



5 Traces of resolvents and spectral densities 

Let us consider quantization of the spinor field on a plane {d = 2) which is orthogonal 
to the Bohm-Aharonov magnetic field configuration. Dirac hamiltonian (16) in gauge 
A — Af, where A/ is given by Eq.(19), takes form 

H = -ia'^dr - ir-^a^{d(p - ie$) + 7°m, (73) 

where $ is the total magnetic flux of the solenoid (Bohm-Aharonov vortex), see 
Eq.(20), and 

a'' = cos if + sin if, a'^ = —a^ sin if + cos ip. 
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In 2 + 1-dimensional space-time the Clifford algebra has two inequivalent irreducible 
representations which can be differed in the following way: 

a^a^^o = is, s = ±1. (74) 

Choosing the 70 matrix in the diagonal form 

7° = c^3, (75) 

one gets 

= -e^''^^'a2e-^''^^% = se^^^^'aie'^^''''^', (76) 

where o"i, o"2, and o"3 are the Pauli matrices, and xi and X-i are the parameters varying 
in interval < Xs < 27r to go over to equivalent representations. Note also that in 
odd-dimensional space-time the m parameter in Eq.(16) can take both positive and 
negative values; a change of sign of m corresponds to going over to the inequivalent 
representation. 

In view of Eq.(74), the kinetic, Eq.(21), and canonical, Eq.(31), angular momenta 
of the planar system take form 

J = -id^ - e$ + ^S7° (77) 

and 

J = -«9^ + ^s7°, (78) 

respectively. 

The kernel of the resolvent (the Green's function) of the Dirac hamiltonian in the 
coordinate representation is defined as 

G-{r,if;r',^') = {r,v\{H - uj)-y,if'), (79) 

where a; is a complex parameter with dimension of energy. Taking into account Eqs.(75) 
and (76), one expands Eq.(79) in modes and gets the following expression 

[r,^,r,^)-^^^2^j l^6„(r;r')e*(^^-^«) c„(r; r')e-(^-^') j ' ^^^^ 

In the case of H given by Eq.(73) radial components of G"^(r, r', yj') ()80|) satisfy 
equations 

/ —uj + m dr + s{n — + s)r^^\ (an{r]r') dn{r;r')\ 

\—dr + s{n — e^)r^^ —uj — m j \bn{r;r') Cn{r;r') J 

_ f —uj + m dr' + s{n — + s)r'^^\ ( an{r;r') bn{r;r') \ _ 

\-dr' + s{n - e$)r'~^ -u - m J \ dn{r; r') Cn{r; r') j 
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All radial components behave asymptotically at large distances (r — oo or r' oo) as 
outgoing waves (e*'^''/ {2TTy/r) or e*'^'" / (27r-\/r')); here k = — and a physical sheet 
for square root is chosen as < ArgA; < tt (Im/c > 0). Another boundary condition 
should be imposed at the boundary of the inner region containing a magnetic vortex. 
In the absence of the vortex, e$ = 0, all radial components obey the condition of 
regularity at small distances (r — * or r' — > 0), and this corresponds to the fact that 
free hamiltonian Hq (i.e. H (73) at 6$ = 0) is essentially self-adjoint. In the presence of 
the vortex, e$ 7^ 0, when a size of the inner region is neglected [r^ 0), the condition 
of regularity at small distances cannot be imposed on all radial components. This is 
due to the fact that hamiltonian H (73) is not essentially self-adjoint, and a parameter 
of the boundary condition at the location of the vortex (at the origin) exhibits itself 
as a parameter of the self-adjoint extension of the hamiltonian operator, for details see 
Refs.pS, 18, 19 . To be more precise, partial hamiltonians are essentially self-adjoint 
for all n with the exception of = ric, where 

nc = m + 1-1^^ (82) 

{uj is the integer part of quantity u (i.e., the largest integer which is less than or equal 
to u). The partial Hamiltonian for n = ric requires a self-adjoint extension according 
to the Weyl - von Neumann theory of self-adjoint operators (see, e.g., Ref . ) . Ap- 
propriately, radial components a„, 6„, c^, and (i„ in Eq.(80) with n ^ are regular 
at r — *• and r' — > 0, whereas those with n = n^. satisfy conditions (for details see 
Ref.lSH): 

cos {^1^ + lirn(|m I r)^a„^(r; r') = — sgn(m) sin {^1^ + lirn(|m|r)^~^6 
cos {^1^ + lini(|m I r)^(i„^(r; r') = — sgn(m) sin {^^1^ + lini(|m|r)"^~^c„^(r; r') 

and 
cos 

cos I 



r; r'^ 



s— H — \ lim (|m|r')'^a„^(r; r') = — sgn(m) sin(s— H — \ lim(|m|r')"^ ^(i„^(r;r') 



e 7r\ ^. , , . . / e 



i^s-^ -|- — j \im {\m\r')^bn^(r;r') = — sgn(m) sin^s-^ + — j lim(|m|r')^ ^Cn^(r;r' 

I 

where G is the self-adjoint extension parameter, and 



B4) 



F = s(e$-Ie$]) + ^-^.; (85) 

note here that Eqs.(83) and (84) imply that < F < 1, since in the case of F = - — -s 

all radial components obey the condition of regularity at r — and r' — 0. Note also 
that Eqs.(83) and (84) are periodic in B with period 2tt. 
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The radial components of the resolvent kernel in the presence and in the absence 
of the vortex are listed in Appendix C. 
Let us consider quantities 



trG'^{r,(p;r',(p) =— ^ [anir;r') + Cn{r;r')] , 

n=-oo 

trJG'^{r,ip;r',<f)=— ^ f n - e$ + -sj [a„(r; r') + c„(r; r') 



n=— oo 
oo 



1 / 1 

trJ'^G'^{r,ip;r',(p)=— J2 [n-e^ + -s) [anir;r') + Cnir;r') 



(86) 
(87) 
(88) 



In the absence of the vortex, using Eqs.(C.14)-(C.17) and performing summation over 
n, one gets in the case of Imk > \Iiek\: 



trC^ir,^; r' = -Ko{-tk\r - r'\), 
tr JoG^(r,(^; r',^) = —Ko{-zk\r-r'\), 



K{){—ik\r — r'\) — — -Ki{~ik\r — r'\ 



(89) 
(90) 
(91) 



where Kp{u) is the Macdonald function of order p. In the presence of the vortex, quan- 
tities (86)- (88) consist of two pieces: one is finite in the limit r' r, and another one, 
which is divergent in this limit, coincides with quantities (89)-(91), correspondingly. 
Moreover, the difference between appropriate quantities in the presence and in the 
absence of the vortex at r' = r is exponentially decreasing as r — >■ oo, and this allows 

2tt oo 

US to perform integration over the two dimensional infinite spatial volume, J d(p J drr, 
and obtain the renormalized traces (see Refs. jHH E2]): 



TrJ{H-Lu)-^ 



^2 _ ^2 



F{u! + m) tan u^i + (1 — F){uj — m)e 



-s{F--)[Tr {H-uo 



tan z/^ + e^P'^ 

sF{l-F) 



,-1 



+ 



^2 _ j^2 





F{l-F)iu 
(92) 

-iF--)u^ + -m 

(93) 

iF(l-F)^-^(F-l)m 



(94) 

here tani/^ is given by Eq.(C.13), and results (92)-(94) are analytically continued from 
region Imk > |Refc| to region Imk > 0, i.e. to the whole complex tu-plane. 

To deal with the divergent at r' — > r pieces of the resolvent kernels, it is sufficient 
to regularize the resolvent kernel in the absence of the vortex and define 



Go'\r,ip; r',ip') = (^r, ip\{Ho - u) ^ exp(-tifo)| r, y?') , 



(95) 
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where t > is the regularization parameter. In Appendix D we obtain the following 
relations in the case of Imk > |Rek|: 



tvGfir,^; r,v) = e~''^^'+''^ E,{-te), (96) 
Ztt 

trJoGfir,^; r, y,) = '-^e-'(-'^'''^E,{-tk% (97) 

An 



(9J 



where 

J u 



du 



is the exponential integral (see, e.g., Ref. [SZj); note that Eqs.(96)-(98) are analytically 
continued to region Imk > 0, i.e. to the whole complex tu-plane. Integrating over the 

2n R 

two dimensional spatial volume, J d(p J drr, where R is the volume radius, we get 



Tr {Ho - c^)-^e-*^o = Ir'u; e-''^' E,[t{m^ - u;% (99) 

TrJo {Ho - cu)-ie-*^o = ^Rhm e-''^' E^[t{m^ - uj% (100) 

Tr Jl (i/o-^)-'e-*^o = \r^uj {RH'^e''"'" + [R^oo' - m") + l]e-'^" E^[t{m' - uj^)]] . 

8 

(101) 

Although traces (99)-(101) are divergent in the limit t 0+, the divergences do not 
contribute to the physical quantities, e.g., Eqs.(63) and (66). This is due to a specific 
form of a discontinuity of the exponential integral at negative real values of its argu- 
ment, ImEi{—u =F ^0) = iivr (m > 0). Consequently, we get finite spectral densities 

Tf\E) = ± lim -ImTT{Ho -Et iO)-^e-*^o = -R^\E\ 9{E^ - m^), (102) 

*->0+ TT 2 

Tf\E) = ± lim -ImTtJo{Ho- ETiOy^e-^^o = -R^ smsgii{E) e{E'^ - m^), (103) 

rf^E) = ± lim -ImTiJl{Ho-ETiQy^e~'^" = -R^\E\[R^{E^ -m') + l]e{E^ -m"). 

(104) 



6 Thermal averages, fluctuations, and correlations 

As follows from two preceding sections, the expressions for thermal characteristics 
consist of two pieces, see Eqs.(61) and (62): the one, denoted by superscript '^^^ and 
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corresponding to the ideal gas contribution, depends essentially on the size of the 
system (R), see Eqs.(102)-(104), and the other one, denoted by superscript and cor- 
responding to the correction due to interaction with a magnetic vortex, is independent 
of R in the limit R ^ oo, see Eqs.(92)-(94). However, it appears that the ideal gas 
contribution to some characteristics is vanishing, and these characteristics are finite as 
the size of the system increases. 

In particular, this is the case for the average fermion number of the system. Let 
us denote the operator of fermion number by N, then it is given by Eq.(40) with N, 
I and 1 substituted for M, J and j\, respectively. Taking into account Eq.(102), one 
gets A^(°)(T) = 0, and using Eq.(92), one gets (see Ref.|31|): 



N{T) 



sin(F7r) 



du 1 /I /I / 

— , tanh -pmyu + 1 



Fu^A - (1 - F)u^-^A-^ + M [(f - I) (u^A + u^-^A-^) - cos(F7r) 

X z — =— : — rr— — : -— — ttt: — — — : — o , „ — ; 



[u^A - ni-^A-i + 2 cos(F7r)]2 + 4(m + 1) sin^(F7r) 



-^9i- cos 9) tanh Q/^E^s ) , (105) 



where 



A 



,1-2F 



rri 



TT 



rx^'"H'2 + 4 



(106) 



T{u) is the Euler gamma function, Ebs is the energy of the bound state in the one- 
particle spectrum, which is determined as a real root of algebraic equation (for details 
see Ref.^) 

note that the bound state exists at cos < (A < 0), and its energy is zero aX A = —1, 
and, otherwise, one has < \Ebs\ < |^| and 



sgn{EBs) = ^sgn(m)[sgn(l + A ^) - sgn(l + A)] . 



Contrary to the average, the quadratic fluctuation of fermion number consists of two 
contributions. Using Eq.(102), one gets 



AW(T;iV,iV) = ^ 



/3|m| 



+ ln(l + e 



-f3\m\ 



(109) 



and, using Eq.(92), one gets 
AW(T; N, N) 



sin(F7r) du , n 
^ ^ ' — sech"' 



271 Jo u \2 
Fu^A + (1 - F)u^~^A~^ - u{2F - 1) cos(F7r) 
[u^A - u^-FA-^ + 2 cos(F7r)]2 + 4(m + 1) sin2(F7r) 



+ 



+ l^(_cose) sech^ 



- -F(l 
4 ^ 



F) sech^ ( -/3m 



(110) 



Let us turn now to another conserved observable - angular momentum. The opera- 
tor of this observable is given by Eq.(40) and, in this instance, we take kinetic angular 
momentum (77) as operator J in Eq.(40). Using Eq.(103), one gets 

jVf(o)(2^) ^ + g-/3H)^ (111) 

and, using Eq.(93), one gets 

mW(T) = -s{F - ]^)N{T) + - F) tanh(^/3m), (112) 

where N{T) is the average fermion number, see Eq.(105). Contrary to the average 
angular momentum, the correlation of angular momentum with fermion number is 
finite in the infinite volume limit: 

A(T; M, N) = ~s{F - i)AW(T; N,N)-^{F- ^-)F{1 - F) sech'{^(3m), (113) 

where A(i)(T; N, N) is the finite piece of the fermion number fluctuation, see Eq.(llO); 
note that the latter relation follows from Eq.(93) and equality A^^\T; M, N) = 
following in its turn from Eq.(103). The quadratic fluctuation of angular momentum 
increases as the volume squared in the large volume limit: 



AW(T;M,M) = ^ 



op2 °? -1 

J ciMMln(l + e-") + (i?W + -)ln(l + e-^H) + 



/3\m 



1 3\m\ 



2e/3|'"l + 1 

where Eq.(104) is used. Using Eq.(94), one gets finite piece 



;ii4) 



AW(T;M,M) = {F--fA^^\T;N,N)--F''{l~Ffsech\-f3m). (115) 

2 8 2 

Turning to nonconserved observables, see Eq.(52), one can take T to be proportional 
either to x ■ o; or to x x a. One can verify that the spectral density in the first case 
is identically zero, whereas it is nonzero in the second case. To be more precise, let us 
take 

T = — X X CK, (116) 
An 

then operator 0(52) corresponds to the observable with the physical meaning of the 
induced magnetic flux multiplied by e and divided by 27r (for details see, e.g., Refs. 
1^). The thermal average of this observable is finite in the infinite volume limit 

e'^sF{l - F) [sin(F7r) J du 
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OCT) = ^ ^ <^ — ^ ^ / — ; tanh -pmy/u + 1 x 

Arcm n J u^u + 1 ^ ^ ' 



X 



u 



l-F A-l 



A-^ - 2ncos(F 



[u^A - u^-^A-^ + 2 cos(F7r)]2 + 4(u + 1) sm^{F 



+ 



+^(-cose) 



m 



(2F-l)EBs + m 



1 



tanh -j3EBs ) + - [F - - ] tanh -/3m 



(117) 



as well as its correlation with fermion number 



A(T;d,iV) 



e^sF{l -F) \ sin(F7r) Z""^ du 



sech 



^l3my/u + 1^ X 



X 



[u^A - u^-^A-^ + 2 cos(F7r)]2 + A{u + 1) sin2(F 



0(-cose)- 



m 



{2F-l)EBs + m 
and its correlation with angular momentum 

A(r;6,M) = 



sech^ (^iiEs^ - \ sech^ (^/3m) | (118) 



A(T;6,iV), 



;il9) 



see Ref.j 

Other nonconserved observables are two pieces of angular momentum, i.e. spin and 
orbital angular momentum. The ideal gas contribution to average spin coincides with 
Eq.(lll), thus the average orbital angular momentum is finite. Expressions for the 
averages of spin and orbital angular momentum and their correlations with conserved 
observables are given in Ref . |32j . 

It should be emphasized that the average of induced fiux and its correlations with 
conserved observables are finite in the infinite volume limit. The average of fermion 
number and its correlation with angular momentum are also finite in this limit, whereas 
the fiuctuation of fermion number, as well as the average of angular momentum, di- 
verges as i?^, and the fluctuation of angular momentum diverges as R^. Thus the ideal 
gas contribution to the thermal characteristics of angular momentum is predominant, 
unless temperature is zero. 

In the limit T — > (/? — * oo) the average of angular momentum tends to finite 
value: 



M(0) 





^n 


'm)(l — 






-1 < A < oo 




iSS^ 


^n 








=0,-1 




iSS^ 


^n 


»[2F2 


-(1- 




-oo < A<-\ 






rn 








-1 < A-^ < oo 


i^s^ 


rn 


^m)(l — 


F)\ 




A = 0,-1 




rn 


»[2(1 


-Ff 




-oo < A-^ < -1 



0<F<i, 



\<F <\, 



(120) 
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whereas its fluctuation tends to zero for almost all values of O with exception of the 
one corresponding to the zero bound state energy, Ebs = (A = —1): 



A(0;M,M) 



0, A^-1, 
\{F-\f, A = -l- 



;i2i) 



the behaviour of the fermion number fluctuation and the correlations in this limit is 
similar to Eq.(121), differing in the values a,t A = —1. It is instructive to present the 
average angular momentum in the form 



M(T) = M(°)(T) + M(0) + m'i^{T), 



'122) 



where M^^\T) and M(0) are given by Eqs.(lll) and (120), respectively, and 



'(1) 



COS 6)- 



sgn(E 



BS) 



exp( 



1 /Stti f 1 
-s{F — -)— — / (it(;sech^(-/9mw) arctan 

2 An J 2 

1 



Ebs\) + 1 
'A{w^- 



1 

r 



sgn(mj 



-A-\ 



exp(/9|m|) + l 
w2-l)i-^ + 2cos(F7r) 



2w sin(F7r) 



(123) 



Thus, averages, as well as fluctuations and correlations, tend exponentially to their 
finite zero-temperature limiting values. 

In the high-temperature limit the averages of fermion number and induced flux tend 
to zero by power law, see Refs. IS21 respectively, whereas correlations tend to flnite 
values: 



A(oo;M,iV) 



-3|(F-i)(l-F)(3-2F), A-VO 
-UF-\)F{l + 2F), A-i = 

-3|(F-i)F(l + 2F), A^O 



< F < i 



(F-i)(l-F)(3-2F), A = 



2' 



(124) 



A(oo;d,iV) = ^'^^^^-^U ^(-cos9), 1 + 

^ ' ' ^ Svrm \ ^ ' {2F -l)EBs + m 2 

sin(F7r) Jdu u^A + u^~^A~^ 1 

^ ^ i V[M^A-ni-^A-i + 2cos(F7r)]V4(u + l)sin2(F7r) J ' ^ ^ 

A(oo; O, M) is proportional to Eq.(125), see Eq.(119). Meanwhile, the average angular 
momentum, the fermion number fluctuation, and the angular momentum fluctuation 
increase as T, T^, and T^, in this limit; such a behaviour is obviously due to the ideal 
gas contribution, see Eqs.(lll), (109), and (114). 

Concluding this section, let us take canonical angular momentum (78) in the ca- 
pacity of total angular momentum of the planar system, then the appropriate operator 
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in the second-quantized theory equals to M + e$iV. It is evident, how to obtain the 
thermal characteristics of this observable from our previous results. In particular, the 
increasing at large R piece of the average remains unchanged, Eq.(lll), whereas the 
finite piece of the average takes form 

M(i)(T) + e$A^(T) = (|e<I)] + ^)iV(T) + - F) tanh(^/5m); (126) 

similarly, the correlation with the induced flux takes form 

A(T; 6, M + e$iV) = ([e$] + ^)A(T; O, N). (127) 

The correlation with fermion number consists of two pieces: 

A(°) (T; M + e$iV, N) = e^A^^^ (T; N, N) (128) 

and 

AW(T; M + e$iV, N) = ([e$] + ^)A«(T; iV, N) - ^{F - i)F(l - F)sech2(l/3m), 

(129) 

as well as the quadratic fluctuation does: 

A(°)(T; M + e$iV, M + e$iV) = (e$)2A(°)(T; iV, N) + A^°\T; M, M) (130) 

and 

A(i)(T; M + e$iV; M + e$iV) = (|e$] + ^)A(^)(T; iV, iV) - 



-FH-F) 



-(F--)e^ + -F(l-F) 



sech2(l/5m). (131) 



It is straightforward to obtain zero-temperature and high-temperature limits of the 
above relations; in particular, the zero-temperature limit of Eq.(126) was first obtained 
in Ref-EU. 



7 Discussion 

In the present paper we consider fractionalization of angular momentum around a 
magnetic vortex in the framework of quantum field theory at finite temperature; for 
the sake of completeness, the results for fermion number and induced magnetic flux 
are also included. If the kinetic definition of angular momentum is chosen, then all 
thermal averages, fluctuations, and correlations are periodic in the value of the vortex 
flux, i.e. they depend on the fractional part of 6$. If the canonical definition of angular 
momentum is chosen, then those thermal characteristics which involve it depend on 
both the fractional and integer parts of e<l>. The difference between two definitions 
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becomes especially striking in the case of the correlation of angular momentum with 
fermion number: the kinetic definition yields finite quantity, see Eq.(113), whereas 
the canonical definition yields a piece which increases in the large volume limit and is 
proportional to e$, see Eq.(128); the latter looks rather unnatural, since such a piece 
should correspond to the ideal gas contribution. The same unnaturalness is relevant 
for the fluctuation of the canonically defined angular momentum: its increasing with 
the volume size piece which should correspond to the ideal gas contribution contains a 
term which is proportional to (e$)^, see Eq.(130). 

According to the generally accepted paradigm, physical manifestations of the Bohm- 
Aharonov effect depend exclusively on the fractional part of e$ (see, e.g., Refs. 
and this favours the kinetic definition of angular momentum. Moreover, the canonical 
definition of angular momentum yields rather embarrassing results for its fluctuation 
and correlation with fermion number. 

Our analysis has been carried out for the whole variety of boundary conditions at 
the location of the vortex, and these conditions are specified by self-adjoint extension 
parameter 0. The nonvanishing of the angular momentum fluctuation signifies that 
angular momentum is not a sharp quantum observable and has to be understood as a 
thermal average only. In the high-temperature limit, angular momentum increases as T, 
see Eq.(lll), and its fluctuation increases as T^, see Eq.(114). In the zero temperature 
limit, angular momentum becomes a sharp quantum observable with finite value M(0) 
(120). However, the last statement is true for all values of in the case of F = i, 
whereas in the case of F 7^ | it is true for almost all values of with the exception 
of one corresponding to the zero bound state energy, Ebs = {A = —1), since in the 
latter case the zero-temperature fluctuation is nonzero, see Eq.(121). 

Among the whole variety of boundary conditions, let us choose the condition of 
minimal irregularity, i.e. the condition corresponding to the radial components being 

divergent at r — at most as r' -p withp < I [laiiHiiini: 

{sf(mod27r), 0<F<| 
0(mod27r), F = \ , (132) 

-sf(mod27r), i<F<l 

or A-^ = OatO<F<i, A = latF = |, A = Oat|<F<l. Under this condition 
thermal characteristics take rather simple form: 

M^'\T) = - 1^ - ^|)'tanh(l/3m), (133) 



A«(T;M,M) = -l(l-|F-i 



2r 1 



F--\+F{l-F) 



sech2(i/3m), (134) 



A(T; M, N) = -i. (f - 1) Q - |F - i|) (1 - |F - i|) sech^(i/3m), (135) 

A{T;d,M)=0. (136) 
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We recall that the results of the present paper are relevant for planar fermions in 
an irreducible 2x2 representation of the Chfford algebra in 2 + 1-dimensional space- 
time, see Eqs.(75) and (76); the mass of such fermions violates parity. As it should be 
expected, our results remain invariant under transitions to equivalent representations 
of the Clifford algebra, i.e. are independent of Xs- As to a transition to the inequivalent 
representation (s — > — s or m — > — m), the averages of fermion number and angular 
momentum are odd and the average induced flux is even under this transition, see 
Eqs.(105), (111), (112), and (117); the quadratic fluctuations of fermion number and 
angular momentum, as well as their correlation, are even and the correlations of induced 
flux with fermion number and angular momentum are odd under this transition, see 
Eqs.(109), (110), (113)-(115), (118), and (119). 

Finally, let us discuss consequences for planar fermions with parity-conserving mass; 
such fermions are assigned to a reducible 4x4 representation which is composed as a 
direct sum of two inequivalent irreducible representations. Thermal characteristics in 
the reducible representation are obtained by summing those in the irreducible one over 
s = ±1. In particular, the average induced flux takes form 



04x4(T) = 2 0(r) 



F=e<I>-[e*] ' 

s=l 

m=\m\ 



(137) 



where 0{T) is given by Eq.(117). The averages of fermion number and angular mo- 
mentum, as well as their correlations with induced flux, are vanishing, 

A^4x4(T) = M4x4(T) = 0, A4x4(T; 6, N) = A4x4(T; 6, M) = 0, (138) 

while the fluctuations of fermion number and angular momentum, as well as their 
correlation, are nonvanishing. 



A4x4(r; 7V,7V) = 2A(T; TV, TV) 



F=e#-[e<l>l 

m—\'ni\ 



(139) 



A4x4(r; M,M) = 2A(T; M, M) 



A4x4(T; M,7V) = 2A(r; M, TV) 



s=l 



= |m| 



s=l 

m=\m\ 



(140) 



(141) 



Note that Eqs.(137) and (141) are flnite in the large volume limit, thus depending 
essentially on the fractional part of e$, whereas the dependence of Eqs.(139) and (140) 
on the vortex flux is not essential, since the ideal gas contribution to fluctuations is 
nonzero and prevailing in this limit. Obviously, our arguments in favour of the kinetic 
deflnition of angular momentum remain valid in this case also. 
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Appendix A. Canonical and improved definitions of 
the moments of momentum tensor of classical fields 

Let T^*^ be the canonical tensor of energy and momentum of classical fields, and let us 
define tensor 

T'*'^ = f ^'^ + aAx''"^, (A.l) 

where x'^'^^ is an arbitrary third-rank tensor which is antisymmetric in two last indices, 
= -X^^". Then 

d^T"" = (A.2) 

and conservation of T^^ leads to that of T'^'^. Moreover, both tensors yield the same 
energy-momentum vector 

P"^ J d^xfi"^^ j d^xT^'^, (A.3) 

since the contribution of transformed into an integral over a surface enclosing 

space, 

J d'x dxx"'' = / d'x dix^'' = j daY'^, (A.4) 

and the latter integral vanishes under mild assumptions on the decrease of x'^^'' at large 
distances. 

Similarly, taking the canonical tensor of densities of moments of momentum 

M'^^' " = x^'f"" - x^'f^" + S"''' (A.5) 

one defines tensor 

^^v-p ^ ^„v;p ^ Q^i^^i^^^P^ _ x^X^"^), (A.6) 

which satisfies divergence relation 

dpM""'" ^dpM^"''^, (A.7) 

and conservation of M'^"''' leads to that of M'*'"'''. 

In cases when tensors T^^'^ and M^'^'^ fail to meet some plausible requirements, the 
transition to tensors T^^ and M^^'P may allow one to eliminate such failings. In par- 
ticular, the canonical energy-momentum tensor of the electromagnetic field is neither 
symmetric nor gauge invariant. By choosing a concrete form of ■^(j^^^ one can remove 
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both of these shortcomings. Incidentally, a convenient expression which does not in- 
volve spin part S^'^'P is obtained for the moments of momentum tensor (Belinfante's 
theorem |2H1)- 

Namely, first, Eq.(A.6) with the use of Eqs.(A.l) and (A. 5) is recast into the form 
M"''' P = x^T^P - x^T^^P + x""^ - X'""" + S''"' P. (A.8) 
Then defining x-tensor as 

^t^uX ^}_^_gl,u,X ^ gu\■,^, _ gXt,;u^ (A.9) 

(recall that S-tensor is antisymmetric in the first two indices and x-tensor is antisym- 
metric in the last two indices), one gets 

MP^'P = xPT^P - x^TPP. (A.IO) 

Thus one gets following expressions for the canonical tensor of moments of momentum 

MP" = J £x [xPT""^ - x^TP^ - dxixPx''^^ - x'^x"^^)] (A.ll) 
and the improved one 

MP" = J d^x {xPT"^ - x"TP^). (A. 12) 

Although the difference between Eqs.(A.12) and (A.ll) is transformed as previously, 
see Eq.(A.4), into an integral over a closed surface, 

j d^xdx{xPx"^^ - x^x""^^) = jda\xPx''''^ - a;"x''°'), (A.13) 

this integral may appear to be finite for certain, long-range, field configurations, since 
it contains an additional power of large distance, |x|. In this case one has to decide, 
whether Eq.(A.ll) or Eq.(A.12) gives physically meaningful moments of momentum. 
Let us consider a system characterized by lagrangian 

C = ~FP"F,,-APj,, (A.14) 

where F^y = d^Ay — d^A^ is the electromagnetic field strength, A^ is the appropriate 
vector potential, is the charge current. The improved energy-momentum tensor of 
the system is proved to be [HH^ (see also, e.g., Ref. |34j) 

TP" = i gp"F^PFxp + FP^Fx" + gP^A^jx - APf (A. 15) 

(note that purely electromagnetic part is symmetric and gauge invariant); incidentally, 
one gets 
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Defining an axial vector from spatial tensor components 

M' = ]^e'^'Mku (A.17) 

and using Eq.(A.12), one gets the following expression for the improved angular mo- 
mentum 

M = j d^x [(x X (E X B)) - (x X A)/], (A.18) 
where E'^ — and — —\e^^^Fki are the electric and magnetic field strengths. 

Appendix B. Derivation of representation (45) for 
thermodynamic potential 

The partition function of the fermionic system is presented as the functional integral 
over the Grassman fields 

exp[-/31](/?,/x)] = j di/j^di/je-^, (B.l) 

where 

Jdrjd'^x i^+{dr -nJ + H)i^ (B.2) 



is the Euclidean action, r is the imaginary time. The integral in Eq.(B.l) is of the 
Gauss type and can be immediately computed 

exp[-/3Q(/3, li)] = det(a^ -/iJ + H). (B.3) 

Hence the thermodynamic potential is given by expression 

= -hndet{dr-fiJ + H) = 
^--Jdr j d'^xti{x,T\lii{dr - f^J + H)\x,t). (B.4) 

^ 

In the case of a static background, operators H and J are r- independent, and the 
integration over r is performed by using the antiperiodicity boundary condition at the 
ends of the imaginary time interval: 

1 oo 

Q(/3,/x) = -- J2 d''xiY{^\\u{H -liJ (B.5) 

r n=-oo 

where c<;„ = ^(n + |), and summation is over integer values of n. Using the notation 
of functional trace, one gets further 

-1 oo -j oo 

= 'IVln(//-/.J-iu;„) = -^TVlnn[(^-/^^)' + ^'] = 

P n=-oo P n=0 
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/3 



TrlnJ] 



n=0 



LOn J 



+ Trln n 



(B.6) 



n=0 



The second term in the figure brackets is dropped as an irrelevant infinite constant, 
and the infinite product in the first term is computed with the use of relation 



cosh(^)= n [l + a\2n+lY 



n=0 



As a result we get expression (45) for the thermodynamic potential. 



Appendix C. Radial components of G"^(r, r', Lp 

The radial components of the resolvent kernel of H (73) take form (see Ref. |31j). 
type 1 (/ = s{n — ric) > 0): 



On t: r = — u; + m) 



r - r')Hl^Jp{kr)Ji^F{kr') + d{r' - r)Ji^F{kr)H\'Jp{kr' 



(1) 



ITT , 



6„(r; r') = y A; [^(r - r')H['^^_p{kr)Ji_F{kr') + 0(r' - r)Ji+^^F{kr)H\'Jp{kr'] 



r(l) 



(C.l) 



(C.2) 
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Cn(r; r') = y (u;-m) [0(r - r')H'i'^^_p{kr)Ji+i_F{kr') + ^(r' - r)Ji+i_F{kr)H\'^^_p{kr') 

(C.3) 

(C.4) 



r(l) 



ITT , 



dn{r; r') = yA; [e{r - r')HtUkr)Ji+iMkr') + e{r' - r)Ji.F{kr)H\'^^_p{kr') 



r(i) 



type 2 (/' = -s(n - nc) > 0): 



it: 



«n(r; r') = y (cu + m) [^(r - r')H^lp{kr)Ji,+F{kr') + 0(r' - r)Jv+F{kr)H\^lp{kr') 



r(i) 



(C.5) 



6„(r; r') 



27r . 



9{r - r')HP_^^p{kr)Jv+F{kr') + d{r' - r)Jv_^+p{kr)Hl^lp{kr 



(1) ^^.^M 



(C.6) 



Cn(?^;?^') = y(u;-m) [6'(r - r')H\}\^p{kr)Jv^i+F{kr') + 6'(r' - r)Ji, ^iJ^F{kr)H\!'\^p{kr') 



r(l) 



27r 



t^n(r; r') = - yA; [^^(r - r')if?^p(A;r) J;,_i+^(A;r') + e{r' - r)Jv+F{kr)H\^\^p{kr' 



r(l) 



(C.7) 
(C.8) 



type 3 (n = Uc): 

in uj + m 



an,(r;r') 



2 sin z/^ + cos u^e''^'^ 



^6{r — r')H^p{kr)[sini>i^J_F{kr') + cos Jir(/i;r')] + 
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+e{r' - r)[smu^J_F{kr) + cos t/^ Jj.(A;r)]i/i^^(A;r')} , (C.9) 

bnSr;r') = {9{r - r')Hi%{kr)[smu^JMkr') + cosu^Mkr')]+ 

+ e{r' - r)[smu^Ji_F{kr) - cos u^J.i+F{kr)]H^2kkr')} , (C.IO) 

Cnc{r;r') = — — — \e{r - r')H^llp{kr){sviiv^J^_p{kr') - co^v^J _x^p{kr')\-^ 

I sm + cos v^e^'^^ ^ 

+ Q{r' - r)[sini/^Ji_F(A;r) - cosi/^ J_i+F(A;r)]//{S(A:r')} , (C.ll) 



— iOir - r')H^l(kr)[smu^Ji_Fikr') - cos u^^ J -i+Fikr')] + 

cosu^e^^'^ 



2 sin + 

+ e{r' - r)[smv^J-F{kr) + cos J F{kr)]H^lF{kr')] . (C.12) 

Here Jp{u) is the Bessel function of order p, H^p\u) is the first-kind Hankel function 
of order p, and 



^^■^ / x/^i ,M-2Fr(l-F) / TT 

" ^7+^'^''^"'''^ '""'^ r(F) ^""^ 2" ^ 4 

In the absence of the vortex the radial components take form: 



(C.13) 



(C.14) 



27r , 



&n(r;r')L^=o = [^(^ - r')HfX,{kr)Jsn{kr') + ^(r' - r) J,„+i(fcr)i/«(A;r') 



(C.15) 

^(r - r')i^i^+i(^r) J,„+i(^r') + e{r' - r)Jsn+i{kT)H^^l,{kr') 

(C.16) 

')"■ 

(C.17) 



Appendix D. Derivation of Eqs.(96)-(98) 



Using integral representation of kernel (95) 



(x (ifo-^)" e-*^o x') 



(27r)^ 



exp 



ip ■ (x — x') — t{p + m ) 



a • p + 7°m + uj 



p2 _ 



we calculate its trace over spinor indices at x' = x 



(x 



2uj 



I 



(Pp e 



-t(p'^ +m?) 



(27r)2 p2_^2 



(D.l) 
(D.2) 
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Similarly, using relations 



Jn (x 



[Ho-uj)-^e-tHl 



exp 



ip ■ (x — x') — t{p'^ + m^^ 



X 



and 



,12 21,1 o\a-p + 7°m + Lo 
X ( xV - xY + -S7°j ^±^, , 



Ho - a;)-^e-*^o | x') = J exp [ip ■ (x - x') - t{p^ + vr?) 



(D.3) 



X 



we find relations 

Jo ( X 

and 



X ) = sm 



(fp e-*(p'+™') 
(27r)2 p2_^2 ' 



(Ho-u;)-'e 



X > = 



d^P r 2 2 



(D.4) 



(D.5) 



(D.6) 



Eqs.(D.2), (D.5), and (D.6) in the case of Imk > \Rek\ are reduced to Eqs.(96)-(98). 
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